The Coefficient of Restitution as a Fluctuating Quantity 
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The coefficient of restitution of a spherical particle in contact with a flat plate is investigated 
as a function of the impact velocity. As an experimental observation we notice non-trivial (non- 
Gaussian) fluctuations of the measured values. For a fixed impact velocity, the probability density 
of the coefficient of restitution, p(s), is formed by two exponential functions (one increasing, one 
decreasing) of different slope. This behavior may be explained by a certain roughness of the particle 
which leads to energy transfer between the linear and rotational degrees of freedom. 

PACS numbers: 45.70.-n, 45.50.Tn, 45.05. +x 



Introduction. The dissipative collision of a solid parti- 
cle with a hard plane may be described by the coefficient 
of normal restitution 



relating the normal components of the relative velocity 
before and after a collision at the point of contact. The 
unit vector n indicates the relative particle positions at 
the instant of the collision. Obviously, s = 1 stands for 
elastic collisions whereas e = indicates the complete 
dissipation of the energy of the relative motion. There are 
several techniques for measuring the coefficient of restitu- 
tion, including high-speed video analysis, e.g., pQ and so- 
phisticated techniques, where the particle is attached to a 
compound pendulum with the axis of rotation very close 
to the center of mass [2 , 3 which makes this method par- 
ticularly suitable for the measurement of the coefficient of 
restitution for very small impact velocities up to mm/sec 
and below. In the presence of gravity, the coefficient of 
restitution can be measured by determining the time lag 
between consecutive impacts of a particle bouncing on a 
hard plane, using a piezoelectric force sensor, e.g., [4j |5] 
or an accelerometer mounted to the plate which detects 
elastic waves excited by the impact, e.g., [6]. When both 
particle and plate are metallic, the time of the impacts 
can be determined by applying a DC voltage between 
ball and plate and determining the instant when the cir- 
cuit closes [6 . In many papers, the time lag and, thus, 
£ is determined by recording the sound emitted from a 
spherical particle bouncing on an underlying flat plane, 
e.g., [7, 8 and many others. From the times U of impacts 
one can compute the coefficient of restitution via 

£{vi) = vj = t^-U^ Vi= 9 {u _ ti _ ih (2) 

where V{ and v[ are the normal pre- and postcollisional 
velocities of the impact taking place at time ti and g is 
the acceleration due to gravity. 

Although it is frequently assumed that the coefficient 
of restitution is a material constant (e.g., this assumption 



is common in many textbooks and widely used in sim- 
ulations), numerous experimental studies have revealed 
that it depends on many parameters: impact velocity, 
material characteristics of the impacted bodies, particle 
size, shape and roughness, and adhesion properties. Here 
we restrict ourselves to the investigation of e of a single 
dry steel ball bouncing on a hard plane, such that the 
coefficient of restitution is only a function of the impact 
velocity e = e(v). 

In several experimental investigations, using either 
photographic techniques [9] or an acoustic emission anal- 
ysis [5], an extraordinary high fluctuation level was no- 
ticed whose origin remained obscure and cannot be at- 
tributed to the imperfections of the experimental setup. 
By means of large-scale experiments as well as microme- 
chanical modeling we make an attempt to characterize 
the fluctuations of the coefficient of restitution and to 
explain the mechanism leading to these fluctuations. 

Experiments. Our experimental setup consists of a 
robot to move a small vertical tube to a desired po- 
sition {x,y,z}. In the beginning of each experimental 
trial a stainless steel bearing ball is suspended at the 
end of the tube by means of a vacuum pump. Switching 
off the pump the sphere is released to bounce repeat- 
edly off the ground where the sound is recorded by a 
piezoelectric sensor. When the ball eventually comes to 
rest, it is pushed to a defined position by a fan where 
it is picked up by the robot who moves the ball again 
to the start position for the next trial. In each cycle 
the initial y}-position is chosen randomly within the 
central region of the ground plate such that edge-effects 
[TP] are not noticeable. The dropping height is chosen 
randomly from z £ [9, 10] cm, corresponding to impact 
velocities v G [1.33, 1.4] m/s which allows for the obser- 
vation of 90-100 bounces of the ball. A massive hard glass 
plate of size 30x20x1.9 cm 3 can either serve directly as 
a ground plate for the experiment or serve as a carrier 
for other ground plates. For our experiments we used 
an extremely hard SiC disk and a Vitreloy disk. The 
measuring process is fully automatized, this allows us 
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to perform a large-scale experiment collecting data from 
thousands of bouncing ball trials. 

From the analysis of the sound-sensor signal we ob- 
tain the impact times ti and, thus, via Eq. ([2| the nor- 
mal impact velocities V{ and the coefficients of restitution 
e(vi). Figure [I] displays the abundance of data, e(y), for 
a stainless steel ball (radius R = 3.0 mm, Young modulus 
Y = 200 GPa, Poisson ratio v = 0.30, density p s = 7.90 
g/cm 3 ). Besides the expected decay of e with increasing 
impact velocity, e.g., [HI [12], we observe a large scatter of 
the data. This scatter becomes apparent only when an- 
alyzing a large number of measurements (here 220,000) 
and was, therefore, not noticed in earlier, similar exper- 
iments [4j [5]. The scatter of s is asymmetric, that is, 



distribution is shaped by a combination of two exponen- 
tial functions (one increasing, one decreasing) of different 
slope. We believe that uncommon statistical properties 
are not due to imperfections of the experimental setup 
but are a consequence of microscopic asperities at the 
sphere's surface. This hypothesis is checked by means of 
a numerical simulation. 

Numerical simulations. The procedure of our simu- 
lation is analogue to the experiment; the rigid ball is 
dropped from a certain height h and repeatedly collides 
with a smooth, hard plate. The air drag is neglected and 
the ball is subjected only to gravity. 

The ball is modeled as a composite multi-sphere par- 
ticle (Fig. [2]). The large central sphere of radius R is 
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FIG. 1: (color online). Experimental results. Top: Coefficient 
of restitution e plotted against the normal impact velocity v. 
The data are colored according to the normalized frequency 
of occurrences. Bottom: Histograms of e for impact velocities 
from small intervals, centered around v — 0.3. .1.0 m/s. The 
lines are exponential fits. 

the deviation of the data with e lower than the mean is 
noticeably larger. Figure [I] (bottom) shows the normal- 
ized frequency p(e) of measurements of a certain value 
of e for several small intervals of the impact velocity v. 
Thus, if we consider e as a fluctuating quantity, p(s) may 
be considered as its probability density function. This 
function reveals strongly non-Gaussian behavior, but the 




FIG. 2: Sketch of the particle model and closeup of its surface. 

randomly covered by many (N ~ 10 6 ) tiny spheres of 
different microscopic size, representing asperities. The 
center of each small sphere i of radius Ri is located at 
the surface of the central sphere. Since Ri/R ~ 10 -3 
(mass ratio rrii/m ~ 10 -9 ) we can safely neglect the 
contribution of the small particles for the computation 
of the moment of inertia of the ball. Thus, the ball is 
characterized by its mass m, i?, center of mass velocity 
vq, angular velocity a;, the set of the radii, of the 
asperities and the set of their position vectors, point- 
ing from the center of the large sphere to the asperities. 
When bouncing, the ball may come in contact with the 
plane either through the central sphere or through one 
of the asperities. In both cases the collision is assumed 
instantaneous and inelastic, characterized by coefficients 
of restitution. 

In the intervals between collisions the ball follows a 
ballistic trajectory 

It — t ) 2 

r G (t) = r G (to)+v G (to) (t-t )+ K 2 0) g, t > t (3) 

where to is the time of the preceding collision and g is 
gravity while the angular velocity uj = uoe^ remains con- 
stant. The evolution of a vector p fixed to the particle is 
then given by 

Prot = P cos out + e w (e w • p ) (1 — cos out) + (e w x p ) sin out 
= A(t)P> 

(4) 
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which defines the rotation matrix A [13]. Therefore, the 
particle contacts the ground when the center of the first 
asperity j reaches the height Rj, that is, the time of 
contact t r follows from the condition 



mm 



r G {t c ) + A{t c )rj 



Rj , 



(5) 



where [...] means the vertical component of the argu- 
ment. The vectorial impact velocity at the point of con- 
tact rv is then 



v c = (v c • n) n + (v c ■ t ) t = v G + u x r c 



(6) 



where n and t are the unit vectors in normal and tan- 
gential directions, see Fig. [3] 




FIG. 3: Sketch of a particle collision. For simplicity only the 
impacting asperity (of exaggerated size) is drawn. 

The post-collisional velocity at the contact point is 
given by 

v' c • n = -e (v c • n) , v' c • t = f3 {v c • t ) , (7) 

with the coefficients of normal and tangential restitution 
e and /3. Finally, we compute the post-collisional velocity 
Vq and angular velocity uj' 



Vq - v G = Av G = —P 



m 



uj' — u> = Acj = —rv x P , 
J 



(8) 



where the transferred momentum P is obtained from 
Eqs. Q and 0: 

v' c -v c = Av G + Aoj x r c = — P + - (r c x P) x r c (9) 

777- f / 

with the mass m and the moment of inertia J of the 
particle. 

Using Eqs. ^ and ([sj) we can compute the dynamics 
of the bouncing particle and, in particular, the times and 
velocities of the impacts, while the coefficients of restitu- 
tion, e and P are given. For e as a function of the normal 



impact velocity v c we use the expression for viscoelastic 
spheres [12] 



e(v c ) = l + ^2dAi(v c y 



/io 



(10) 



with Ci being known constants [C\ — Cs = 0; C2 = 
-1.153; C 4 = 0.798; C 5 = 0.267. . . , see for details) 
and Ai being material constants. We determined the first 
three non-trivial constants, A2 ~ 0.0467; A^ w 0.1339; 
A 5 « -0.2876, by fitting Eq. (10) to the experimental 
data shown in Fig. [l] 

The coefficient of tangential restitution j3 depends on 
both bulk material properties and surface properties. 
Therefore, (3 cannot be analytically derived from material 
properties, except for the limiting case of pure Coulomb 
friction [14]. Here we use /3 = 1. 

Simulation results, On the macroscopic level (ne- 
glecting the microscopic asperities at the surface of the 
particle), the simulation results can be analyzed in the 
same way as the experimental data. We introduce the 
macroscopic coefficient of restitution e as the ratio of the 
post-collisional to pre-collisional center of mass velocities 
in normal direction 

e = -^. (11) 

V G 71 

Surprisingly, the macroscopic interpretation of our simu- 
lation results depicted in Fig. [4] shows a striking resem- 
blance with the experimental data. While the actual co- 
efficient of normal restitution e is a function of the normal 



impact velocity described by Eq. (10), the macroscopic 
coefficient of restitution e computed via Eq. ( [TT] ) reveals 
strong fluctuations in agreement with the experiment. 

The probability densities, p(e), (bottom panel of Fig. 
|4| are also close to the distribution obtained in the ex- 
periment. Their shape is excellently approximated by a 
combination of two exponential functions, one increas- 
ing and one decreasing. The peaks of p(e) are in line 
with the values of s(v G ) obtained from Eq. (10) for the 
corresponding velocities. 

For the case of an absolutely smooth particle (no as- 
perities), the impact velocity v c at the point of contact 
is the same as v G and e is equal to the coefficient of nor- 
mal restitution e. Therefore, in this case, there would be 
no scatter of the measured values of e. Consequently a 
scatter of these data must be attributed to the asperities 
at the surface of the particle. From the agreement of the 
experimental data, Fig. [l] and the simulation data, Fig. 
[4j it became apparent that the tiny microscopic imper- 
fections of the surface of the otherwise macroscopically 
smooth particle lead to the characteristic fluctuations of 
the coefficient of normal restitution. The slight asymme- 
try in the shape of the probability density in case of the 
experimental results, i.e., p(s) compared to may be 
caused by some additional dissipative forces that, how- 
ever, have not be taken into account in our model. 
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FIG. 4: (color online). Simulation results. Top: The macro- 
scopic coefficient of restitution e plotted against the normal 
center of mass velocity vg in moment of impact. The data are 
colored according to the normalized frequency of occurrences. 
Bottom: Histograms of s for velocities from small intervals, 
centered around i;g = 0.3..1.0m/s. The lines are exponential 
fits. 



Let us consider the role of the sphere's rotational de- 
grees of freedom which may be interpreted as internal 
degrees of freedom since they do not enter the computa- 
tion of the coefficient of restitution, neither in the exper- 
iment, Eq. (J2|, nor in the simulation, Eq. (11). Initially, 
the particle had no spin and only one degree of freedom 
in its translational motion. However, due to eccentric 
impacts caused by asperities, the particle gains some ro- 
tation and acquires velocity in horizontal direction. The 
partition of total kinetic energy in rotational and trans- 
lational degrees of freedom constantly varies from one 
impact to another. Thus, the kinetic energy of the lin- 
ear vertical motion (the only component which enters e) 
just before a collision is transformed into energy of the 
rebound vertical velocity, dissipated energy due to the 
coefficient of restitution, and changes of the horizontal 
and rotational velocities. The latter two contribution 
may be positive or negative, leading to a reduced or in- 
creased values of the measured coefficients of restitution, 
according to Eqs. Q and (11). Therefore, the particle 



rotation may be considered as a reservoir of internal en- 
ergy, leading to fluctuations of the measured coefficient 
of restitution. 

The substantial increase of data scattering with the de- 
crease of impact velocity can be attributed to the grow- 
ing role of the rotational degrees of freedom in the energy 
partition. In consecutive collisions the translational en- 
ergy decreases appreciably due to dissipation, but the 
amount of energy concentrated in rotation varies only 
slightly. Thus, the proportion of rotational energy to 
translational energy increases with the number of parti- 
cle bounces. The kinetic energy transfers from rotational 
to translational mode and vice- versa result in a more ap- 
parent fluctuation of the coefficient of restitution. As a 
consequence, the scattering of the restitution coefficient 
increases as the impact velocity decreases. 

Conclusion. We have performed an experimental and 
numerical study of the coefficient of normal restitution 
as a function of impact velocity. From about 2.2 x 10 5 
experiments of the same stainless steel sphere bouncing 
on a massive horizontal plate we determined experimen- 
tally the probability distribution of the fluctuations of 
the coefficient of normal restitution, e. We found that s 
increases as the impact velocity decreases. For fixed im- 
pact velocity, the probability density of the coefficient of 
restitution p(e) is non-Gaussian. It consists of two expo- 
nential functions (one increasing, one decreasing) of dif- 
ferent slope. We modelled the particle used in the exper- 
iment by a mathematical sphere whose surface is covered 
by a large number of much smaller spheres (asperities) to 
simulate a certain roughness. The simulations revealed 
the same properties of the fluctuations. Since the as- 
perities are the only origin of scatter, we conclude that 
the experimental observed fluctuations of the coefficient 
of restitution coefficient are due to microscopic surface 
roughness of the ball, causing energy transfer between 
the translational and rotational degrees of freedom. 
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lence Engineering of Advanced Materials. 
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